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Abstract 



We examine two different ways of encoding a counting function, 
as a rational generating function and explicitly as a function (denned 
piecewise using the greatest integer function). We prove that, if the 
degree and number of input variables of the (quasi-polynomial) func- 
tion are fixed, there is a polynomial time algorithm which converts 
between the two representations. Examples of such counting functions 
include Ehrhart quasi-polynomials, vector partition functions, integer 
points in parametric polytopes, and projections of the integer points in 
parametric polytopes. For this last example, this algorithm provides 
the first known way to compute the explicit function in polynomial 
time. We rely heavily on results of Barvinok, and also of Verdoolaege, 
Seghir, Beyls, et al. 



1 Introduction 

We are interested in a wide variety of functions of the form 



Most examples, including Ehrhart quasi-polynomials and vector partition 
functions, will count some combinatorial object. The function c(s) can be 
encoded in at least two different ways: either as an explicit function or as a 
generating function 



c : Z 



•n 



Q 
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Example 1.1. Consider the generating function 

/(x) = -^ = l + x 2 + a; 4 + --- = ^ 



c(s)x 



The corresponding function can be represented explicitly as 

0, if s < 
c(s) — { 0, if s > and s odd 

1, if s > and s even. 



□ 



Mathematicians often encode a function as a rational generating function, 
such as f(x) = in Example II. 1[ which is a compact representation of a 
(possibly infinite) Laurent power series J2 s ez n c ( s ) xS > where c(s) e Z n . This 
has the advantage that we may apply many computational tools to manipu- 
late our rational generating function and obtain information from it (see, for 
example, |BW03j ). An explicit function representation for c(s), on the other 
hand, has the advantage of being easily evaluated for a particular value of s. 
Such a representation is therefore preferred in the compiler community (see, 
for example, [VSB + 07j ). 

We will show that these ways of representing a function are "the same," 
in the sense that one can convert between the rational function and explicit 
function representations in polynomial time (if the degree and number of 
variables of the function is fixed). Let us be more precise about the specific 
representations we will use for generating functions and explicit functions. 

Definition 1.2. By a rational generating function /(x), we will mean a 
function given to us in the form 



/(x) = V a t — r- — , (1.3) 

(l-x b ")(l-x b ^)---(l-x b ^) 1 ; 



OLj 

Vl - vbiiVl - 

i6/ 

where x e C™, J is a finite set, e Q, G Z n , and bjj {0}. 

Definition 1.4. A step-polynomial g : Z n —>■ Q is a function written in the 
form 

m dj 

9( s ) = Yl a i II L( a i^ s ) + h ik\ , 
j=i k=i 
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where aj G Q, a jfc G Q n , bj k G Q, (•, •) is the standard inner product, and |_-J 
is the greatest integer function. We say that the degree of g(s) is max.,{ dj }. 

A piecewise step-polynomial c : Z n — ► Q is a collection of polyhedra Qi 
(which may not all be full dimensional) together with corresponding functions 
g t : Qi n Z n -> Q such that 

1. the int(Qj) partition Q n (where mt(Q) is the relative interior of Q in 
the affine space it lies in) 

2. c(s) = flfj(s), for s G int(Qi) H Z n , and 

3. each g { is a step-polynomial. 

We say that the degree of c(s) is maxjdeg^j. Working with the relative 
interiors of the polyhedra allows us not to worry about the value of the 
function at the intersection of two polyhedra. 

For example, the explicit function c(s) in Example 11.11 can be written as 
the piecewise step-polynomial 

fi + LiJ-L^J> if *>o 

c(s) = < 1, if s=0 

[o, if s < 0. 

We must be careful when speaking of a correspondence between a rational 
generating function and a piecewise step-polynomial, because a generating 
function may have different Laurent power series expansions which converge 
on different regions of C n . For example, if f(x) = then 

1 + x + x 2 + x 3 + ■ ■ ■ and — x~ x — x~ 2 — x~ 3 — ■ ■ ■ 

are Laurent power series expansions convergent for ||x|| < 1 and > 1, 
respectively. 

We state the main theorem, and then provide several examples of rational 
generating functions and piecewise step-polynomials. 

Theorem 1.5. Fix n and k. There is a polynomial time algorithm which, 
given a rational generating function /(x) in the form U.3\) with n variables 
and each hi < k and given 1 G U 1 such that (1, b^) 7^ for all i and j, 
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computes the piecewise step-polynomial c : Z n — > Q tyii/i degree at most k 
such that 

/(x) = £ c(s)x s 

is t/ie Laurent power series expansion of /(x) convergent on a neighborhood 
of e 1 = (e' 1 , e' 2 , . . . , e'") ; e £/ie frase 0/ i/ie natural logarithmic function. 

Conversely, there is a polynomial time algorithm which, given a piece- 
wise step-polynomial c : Z n — > Q of degree at most k such that /(x) = 
SseZ" c ( s ) xS converges on some nonempty open subset ofC n , computes the 
rational generating function /(x) in the form U.3\) with ki < k. 



The proof of the first half of this theorem will use several ideas from 
[BP99]. Section 3 will be devoted to the proof of the theorem, after we lay 
the groundwork in Section 2. Note that applying the theorem twice (in one 
direction and then the other) will in general not result in the exact same rep- 
resentation of the rational generating function or piecewise step-polynomial. 
We are unaware of any canonical form for either rational generating functions 
or piecewise step-polynomials that can be computed in polynomial time. 

As there may be many functions with the same generating function repre- 
sentation (convergent on different neighborhoods), we need to find an appro- 
priate 1 value when we want to convert a given rational generating function to 
an explicit representation. If we know that the function c(s) is only nonzero 
for s in some polyhedron Q such that Q does not contain any straight lines, 
then we may take any 1 such that (1, by) 7^ for all i,j and such that 

Q n { x e Q n I (l,x) >o} 

is bounded. Such an 1 will give us the desired Laurent power series expansion 
Y2 S c(s) x s . In Example 11.11 we could take I = —1. 

Example 1.6. Let P C Q d be a rational polytope, and let 

c P (s) = #( s pnz d ), 

where sP is P dilated by a factor of s. 
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Then Ehrhart proved |Ehr62| that cp(s) is a quasi-polynomial, that is, 
there is a V G Z + and polynomial functions go(s), g±(s), . . . , gv-i( s ) such 
that 

cp(s) = gj(s) for s = j (mod Z>). 

The generating function Yl^=o Cp ( s ) xS can ^ e computed in polynomial time, 
and this has been implemented in LattE (see [DLHTY04J). Computing 
some explicit function representation of cp(s) in worst-case exponential time 
has been implemented in |CL98j and computing cp(s) as a piecewise step- 
polynomial in polynomial time had been implemented in |VSB + 07j . 

Example 1.7. In particular, let P C Q 2 be [0, |] x [0, |]. 

Then 



cp(s) 



for s > 0, 



and we have that 



s=0 



(l-x)(l-x 2 ) 2 
which can be verified by hand. 



^1 —x)(l — x 2 ) 



□ 



Example 1.8. Given ai, a 2 , . . . , a^ G N n , let c : Z n — > Z be the vector 
partition function, defined by 

c(s) = # { A = (Ai, A 2 , . . . , A d ) G N d | s = Aiai + A 2 a 2 + ■ ■ • + A d a d } , 

i.e., the number of ways an integer vector s can be written as a nonnegative 
combination of the aj. 

Then the generating function representation of c(s) is very simple: 



v ; (1 - x a i)(l - x a 2) • ■ ■ (1 - x a <*)' 

The piecewise step-polynomial representation of c(s) can also be computed 
in polynomial time (see Corollary 13. ll or |VSB + 07j ). Beck [Bec04j describes a 
general technique for computing vector partition functions, based on partial 
fraction expansions of /(x). He does not provide a complexity analysis, but 
standard techniques for computing partial fractions |Hen74j are exponential, 
even for fixed dimensions. 
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Example 1.9. In particular, consider the number of ways to partition an 
integer s into 2's and 5's, i.e., a\ = 2 and 02 = 5. Then the generating 
function representation is 



/(*) 

and 



1 



;i -x 2 )(i -x 5 ) 



c(s 



0, if s < 

li« + ij + HW > tf«>o, 

which, again, can be verified by hand. □ 



Both Ehrhart quasi-polynomials and vector partition functions are special 
cases of counting integer points in parametric polytopes. In general, we let 
P c Q n x Q d be a rational polyhedron such that, for all s 6 Q", the set 
P s = {t G Q d I (s, t) G P} is bounded, and we define the function c : Z n —>■ Z 
by 

c(s) = #(P S n z d ) = # {t G z d I (s,t) GP}. (1.10) 

We call P a parametric polytope, because, if P = |(s, t) G Q n x Q"' | 
As + Bt < c} for some matrices A G Z mxn , B G Z mxd and vector c G Z m , 
then 

P s = {t G Q d I Bt < c - As} , 

so as s varies, the polytope P s varies by changing the right hand sides of its 
defining inequalities. 

Both a piecewise step-polynomial representation for c(s) and its generat- 
ing function, c(s)x s , can be computed in polynomial time, as the following 
two propositions state. 

Proposition 1.11 ( |VSB + 07j ). Fix n and d. There is a polynomial time 
algorithm which, given a parametric polytope P C Q n x Q d , computes the 
piecewise step-polynomial 

C (s) = #(P S n z d ) 

with degree at most d. 
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Proposition 1.12. Fix n and d. There is a polynomial time algorithm which, 
given a parametric polytope P C Q" x Q d such that 

/(x) = £ c(s)x s 

converges on some nonempty open subset ofC n , computes /(x) as a rational 
generating function of the form U.3\) with the kj at most d. 



In Section [2] we recall the key ideas of the proof of Proposition 11.111 from 
[VSB + 07] . drawing heavily from the ideas in [BP99] . Proposition 11.121 is an 
immediate consequence of [BP99 , Theorem 4.4] and the monomial substitu- 
tion from |BW03] and will be proved in Section [31 

We may also look at projections of the integer points in a parametric 
polytope. Let P C Q n x Q d x Q m be a rational polytope, and define the 
function c : Z n — > Z by 

c(s) = # {t G Z d | 3u G Z m : (s, t, u) G P} . 

If P s = { (t, u) G Q d x Q m | (s, t, u) € P} and the projection vr : Q d x Q m -> 
Q d is defined by 7r(t, u) = t, then 

c (s) = #(7r(P s nz d+m )). 

It follows from [BW03] that the generating function, X] s c ( s ) xS > can ^ e com " 
puted in polynomial time (for fixed n, d, and m). Therefore, we have as a 
corollary to Theorem 1 1 . 51 1 hat the piecewise step-polynomial can be computed 
in polynomial time. 

Corollary 1.13. Let n, d, and m be fixed. There is a polynomial time 
algorithm which, given a polytope P C Q n x Q d x Q m , computes the piecewise 
step-polynomial 

c(s) = # {t G 7L d | 3u G Z m : (s, t, u) G P} 

with degree at most n + d + m. 



We will prove this corollary at the end of Section 3. 
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Computing Piecewise Step- Polynomials for 
Parametric Polytopes 



In this section, we recall the main elements of the proof of Proposition 11.111 
from [BP99] and [VSB+07j . That is, given a parametric polytope P C Q n x 
Q d , we define P s = {t G Q d | (s, t) G P}, for s G Z n , and we want to compute 

c(s) = #(P S n z d ) 

as a piecewise step-polynomial. We demonstrate each step with a running 
example and formulate an extended version of the final step for use in Sec- 
tion [3j 

Example 2.1. Consider the parametric polytope 



P 



(s,t) G 



X 





-1 


2\ 




/ 1 


-2\ 




1 


-1 




-1 


1 










s + 


1 










0/ 









t > 



We want to compute a piecewise step-polynomial representation of 
c(s) = #(P S n Z 2 ) = # {t G Z 2 | (s,t) G P} . 



□ 



Our main tool will be a slightly different sort of generating function than 
we have been using. If S C Z ' is a set of integer vectors, then define its 
generating function to be 

/(5;x) = ^x t = tftf '-fi- 

tes (t!,...,t d )es 

In our previous notation, this is the generating function for c(t) such that 
c(t) = 1 for t G S and c(t) = otherwise. 

Our proof of Proposition 11.111 will have two main steps. 

• First, we will calculate the generating function f(P s fl Z d ; x) as a ratio- 
nal generating function, and we will examine how it changes as s varies 
(Propositions El and ETLTf . 
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• Second, we will calculate 

c(s) = #(P s nz d ) = /(P s nz d ; i), 

by appropriately substituting x = 1. 

In order to calculate the generating function f(P s flZ d ; x), it is necessary 
to know what the vertices of P s are. 

Example 2.2. Consider the parametric polytope P from Example 12.11 

For a given s, the vertices of P s can be obtained as the intersections of 
pairs of facets of P s . The facets t\ — and t\ — 2t 2 = s% — 2s 2 , for example, 
intersect at the point vi = (0, — si/2 + s 2 ). This point is not always active, 
that is, actually a vertex of P s . It is active exactly when 2s 2 > s\ > (for 
all other values of s, vx ^ P s ). We similarly find the vertices v 2 = (0,0), 
v 3 = (si - s 2 , 0), v 4 = (si - 2s 2 , 0), v 5 = (0, -s 1 + s 2 ) and v 6 = (s l5 s 2 ), 
active on the domains 2s 2 > si > s 2 , si > s 2 > 0, si > 2s 2 > 0, s 2 > si > 0, 
and Si,s 2 > 0, respectively. Combining all of the inequalities, we have the 
regions 



Qi 


= {s 


2s 2 > s 1 > s 2 } 


Q 2 


= 


si > 2s 2 > 0} 


Qs 


= {s 


s 2 > si > 0}. 



For s e Qi, the polyhedron P s has active vertices v l5 v 2 , v 3 , v 6 ; for s e Q 2 , it 
has active vertices v 3 , v 6 , v 4 ; and for s e Q3, it has active vertices v 1; v 5 , v 6 . 
On the boundary of the Qi, there is more than one possible description of 
the vertices (any is fine). 

Figure 12.31 shows the decomposition, the vertices active in each Qi, and 
the evolution of the vertices as the value of s changes. 

□ 

As the example suggests, and as shown in |VSB + 07j . we can find poly- 
hedra Qi such that the int(Qi) partition Q n and, for any s in the relative 
interior of a given polyhedron Qi, the polytopes P s will have a fixed set of 
vertices given by affine transformations of s (where an affine transforma- 
tion T : Q n — > Q d is given by T(s) = T'(s) + v such that T' is a linear 
transformation and v G Q d ). These Qi will be the pieces of our piecewise 
step-polynomial. This is the content of the following proposition. 
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Proposition 2.4 (Decomposition). Fix d andn. There exists a polynomial 
time algorithm, which, given a parametric polytope P C Q n x Q d , finds 
polyhedra Qi whose relative interiors partition Q n , and, for each i, computes 
a collection of affine transformations Tn,Ti2, . . . , Ti mi : Q n — > Q d , such that, 
for s G int Q h the vertices of P s are 2a (s), T i2 (s), . . . , T imi (s). 

Algorithms to compute the parametric vertices and the chambers can be 
found in |LW97j and |CL98j respectively. A proof of the polynomial time 
complexity is given in [VSB+071 - 

Now we can concentrate on computing f(P s H Z d ;x), given that s is in 
the relative interior of a particular Qi. As a first step, we examine how 
to compute the generating function of an easy set: the integer points in a 
unimodular cone. The general case of a polyhedron is based on a reduction 
to these unimodular cones. 

Definition 2.5. Let Ci,C2,. . . ,q € Z d be a basis for the lattice Z d , and let 
Pi G Q, for 1 < i < d. We define the rational unimodular cone 

K={xGQ d (a, x) < Pi for 1 < % < d } . 

This cone may have a vertex which is not at the origin. Let u 1; u 2 , . . . , 

be the negative dual basis of Z , so that 



(u,, Cj) 



-1, if i = j 
0, ifz^j. 



If Pi were zero, for all i, then K would be the cone with vertex at the origin 
defined by 

K = { Aiui + A 2 u 2 + ■ ■ • + \ d u d | A>0}, 
we would have that 

K n Z d = { Aiui + A 2 u 2 + ■ ■ ■ + X d u d | A G Z> }, 

and therefore 

f(K n Z d ; x) = — s 7 -. 

JK ' (l-x u i)(l-x u 2)---(l-x u d) 

In the general case, where the Pi are not necessarily zero, we have that 



where p = — ^2 i=x LA J u i ( see [BP99] ). This greatest integer function in the 
definition of p is where the greatest integer function in our step-polynomial 
will come from. Note also that the denominator of this generating function 
does not depend on the A, only on the Cj. 

We want to reduce our problem, which is finding the generating function 
f(P s H Z d ; x) where P s is a polyhedron, to the easy problem of finding the 
generating function for a unimodular cone. We can first reduce to the case 
of (not necessarily unimodular) cones using Brion's Theorem [Bri88] . which 
states that the generating function of a polytope is equal to the sum of the 
generating functions of its vertex cones. These vertex cones are formed by 
the supporting hyperplanes of the polytope that intersect in a given vertex 
(see Figure |2"U1 for an example). Next, we use Barvinok's unimodular de- 
composition |Bar94] to write the generating function of each vertex cone as 
a (signed) sum of generating functions of unimodular cones. 

Example 2.7. Consider once more the parametric polytope P from Ex- 
amples 12.11 and 12 .21 We want to compute the generating function of this 
parametric polytope. Consider specifically region Q% from Example 12.21 with 
active vertices vi = (0, — si/2 + s 2 ), v 5 = (0, — S1 + S2) and v 6 = (s\, S2). The 
polytope corresponding to s = (3, 4) £ C3 is shown in Figure 12.91 together 
with the vertex cones, cone(P s ,Vj), at each active vertex. Brion's theorem 
tells us that 

f{P s nZ 2 ;x) = /( cone(P s , Vl ) fl Z 2 ; x) + /( cone(P s , v 5 ) n Z 2 ; x) 

+ /(cone(P s ,v 6 )nZ 2 ;x). 

The vertex cones at V5 and V6 are unimodular, but the one at vi is not. We 
therefore need to apply Barvinok's unimodular decomposition to cone(P s , Vi), 
which yields 



/(cone(P s , Vl )nZ 2 ;x) 



.(o,-l*-«J) 



(1 - x( 1 '°))(l - x^ 1 )) (1 - x( 1 '°))(l - X^' 1 )) ' 

. ., . (2-8) 

We refer to [BP991 IDLHTY041 |Kop07| IKV07] for details on how to perform 

Barvinok's decomposition. Table 12.101 lists the generating functions of all 

vertex cones. 



□ 
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Vi 








v 6 




Vi 






















Figure 2.9: -P(3,4) and its vertex cones 



Vertex v,- 



/(cone(P s ,v i )nZ d ;x) 



Vl = (0, -8 X I2 + s 2 ) 



t (_2L4J— » 2 J+«i-2«a.-L^— »2J) 
(l-x^Xl-xO 2 . 1 )) 



x (o,-L^- S2 J) 
(l-x( 1 .°))(l-x(°. 1 )) 



v 2 = (0,0) 



c (0,0) 



(l-x( t ». 1 ))(l-x( 1 .°)) 



v 3 = (si - S 2 ,0) 



x ( s l- s 2.°) 
(l-x( 1 .°))(l-x(- 1 .°)) 



v 4 = (s x - 2s 2 ,0) 



( S1 -2s 2 ,0) 



(l-x(2. 1 ))(l-x( 1 .0)) 



v 5 = (0, -si + s 2 ) 



(0,- Sl + S2 ) 



(l-xCWXl-xC 1 . 1 )) 



v 6 = (si,s 2 ) 



C (S1,3 2 ) 



(l-xf- 2 - 1 ))^^- 1 " 1 )) 



Table 2.10: The generating function of each vertex cone 
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When we do this in general, the end result is the following proposition, a 
rephrasing of Theorem 4.4 of [BP99J. 



Proposition 2.11. Fixd. There exists a polynomial time algorithm, which, 
given a parametric polyhedron P C Q n x <Q 
that for s G Q the vertices of P s = {t G Q d 

transformations Ti(s),T2(s), • • ■ ,T m (s), computes the generating function 



d and a polyhedral region Q such 
(s, t) G P} are given by affine 



/(P s nZ d ;x) = ^^ 



r Pi(s) 



1 



rbi] 



)(1 



1 



where e E { —1, 1 } ; bjj G Z d \ {0} ; and eac/i coordinate o/pj(s) 
zs a step-polynomial of degree one, for each i. 



Now that we know how to compute f(P s H Z d ; x), all that remains is to 
evaluate it at x = 1. 

Example 2.12. Consider once more the parametric polytope P from Ex- 
amples I2TT1 [2721 and 12771 

We have already computed f(P s D Z 2 ; x), and we now compute the value 
f(P s fl Z 2 ; 1). We cannot simply plug in x = 1, because 1 is a pole of some 
of the rational functions. Instead, we make a suitable substitution, in this 
case x=(t+l,t + l) (chosen carefully so that none of the denominators 
become identically zero), and take the limit as t approaches zero. To compute 
this limit, we can simply compute, for each term in the sum constituting 
f(P s flZ 2 ;t + l,t + l), the constant term in the Laurent series expansion at 
t = 0. 

For example, substituting x = (t + l,t+l) into the second term in (12. 8p . 
we obtain 

(l + t)-L%]+« 2 
(1 -(! + *))(! -(! + *))■ 

Since the denominator, in this case, is exactly t 2 , the constant term in the 
Laurent expansion is simply the coefficient of t 2 in the expansion of the 
numerator, i.e., 



■LiJ+^x-LiJ+^-i; 



1 


Sl 


2 s 2 




1 

S2 ~2 


Sl 


2 


.2. 


+ 

2 


. 2 . 


. 2 . 



+ 



S-2 



The other terms are handled similarly. Note that this is the place where the 
step polynomials show up in full force. The contribution of each vertex cone 
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Vertex v. 


/( cone(P s , Vj) fl Z d ; l) (if Vi is active) 


Vl = (0,- Sl /2 + s 2 ) 


4 + -*i LfJ-f -f + LfJ 2 +LfJ + I 




~\ Lf J ~ ^ + Lf J 52 ~ \ Lf J + f 


v 2 = (0,0) 





V 3 = (Sl - S2,0) 


2 2 
s l 1 S 1 S 2 s 2 1 1 

4 2 4 8 


v 4 = (si - 2s 2 , 0) 


»l 2sis 2 si | 2s l 1 1 2 
6 3 2~ r 3~ r 2_r 9 


v 5 = (0, -si + s 2 ) 


2 2 
s l «1S2 1 Sl 1 s 2 S 2 _i_ 1 

4 2~ r 2~ r 4 2" r 8 


v 6 = (si, s 2 ) 


2 2 
s l _i_ Sl«2 _i_ «1 _i_ s 2 _i_ S2 i 47 

12 ' 6 ^ 2 ~ r 12 ' 2 ~ r 72 



Table 2.13: The contribution of each vertex cone to the constant term of the 
Laurent expansion 

to the constant term of the Laurent expansion is listed in Table 12.131 The 
final step-polynomial in each chamber is computed using Brion's Theorem 
as the sum of the appropriate step-polynomials from this table. The final 
result is shown in Figure 12. 141 

□ 

In general, we use the following lemma, which is more general than strictly 
needed here, but which allows for an incremental computation as discussed 
after the lemma and which we will also need in Section [3j The lemma is 
a special case of the monomial substitution theorem [BW031 Theorem 2.6]. 
We provide a slightly different proof, which lends itself more easily to an 
implementation. It is an extension of an idea from |DLHTY04j . which is in 
itself a variation of the idea used in [Bar 94] . 

Lemma 2.15 (Specialization). Let us fix k. There exists a polynomial time 
algorithm which, given a rational generating function /(x) of the form U.3\) 
and an m with < m < d such that g(z) := f(zi, . . . , z m , 1, . . . , 1) is an 
analytic function on some nonempty open subset of C m , computes g(z) in 
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\ 




X 



Figure 2.14: The enumerator of P s , a step-polynomial in each chamber 
the same form, i.e., 



where k v <k, z G C m , G Q, q^ G Z m , and d iY G Z m \ {0}. 

Furthermore, if the vectors and i/ie numbers cti in U.3\) are fixed, 
but the vectors pi vary, then the vectors dj/j/ are fixed, each differ by a 
constant vector from some Pi, and fa are each a polynomial of degree at most 
k in the coordinates of some p^. 

Proof. The case m — d is trivial, so we will assume m < d. Note that we 
cannot simply plug in the values 1, since (z±, . . . , z m , 1, . . . , 1) may be a pole 
of some of the terms in (II. 3p . In fact, if m = 0, then it will be a pole of all 
those terms. We must take an appropriate limit as (x m+ i, . . . , x n ) approaches 
(!,...,!). Consider 



as a function of t only, where A G Z d m is such that for each i 6 J either 
(bn, . . . , b im ) ^ or ((6 i>m+1 , . . .,b id ), A) 7^ 0. Such a A can be found in 
polynomial time by choosing an appropriate point from the "moment curve" 




(2.16) 



h(t) = f(z 1 ,...,z m ,(l + t) x \...,(l + t) x *-™) 
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as in |BP99} Algorithm 5.2]. Then g(z) is simply the constant term in the 
Laurent power series expansion of h(t) about t = 0. This is the sum of the 
constant terms in the Laurent power series expansions of 

hi(t) =ai- 



z b »i(t + l)«ii)(l - z b «(t + ... (1 - z b ^ (t + 1) 



: c + cit + c 2 t 2 + 



where, for v e { Pi, by }, we write v' for the first m components of v and v" 
for the remaining d—m components, and we let = (p", 1) and % = (by, 1). 

Consider a particular ^j(t). Let r be the number of factors with fjj ^ 
but = 0. Then hiit) has a pole of order r at t — 0. Therefore, we must 
compute the coefficient of t r in the Taylor series expansion of t r hi(t), which 
is analytic at t — 0. 

Following [DLHTY04] we use the technique outlined in [Hen74t 241-247] 
(where it is applied to compute the residue of a function, i.e., the coefficient 
of the term t^ 1 ). Let t r hi(t) = where P and Q are polynomials. To 
compute the coefficients Cj in 

P(t) 

Qit) 

expand Pit) and Q(t) as 

P{t) =: ao + ait + a 2 t 2 + 
Q{t) =: b + M + ht 2 + ■ 

and apply the recurrence relation 

^ = — I a, - b ' r :, 

\ i=l 

Note that we only need to keep track of the first r + 1 coefficients of P{t) and 
Qit), and so this may be done in polynomial time. Examining the recursive 
process, we see that the lemma follows. □ 

Remark on the implementation of Lemma I2.15t Note that as argued 
by [DLHTY04] , a A from the moment curve may not be the most appropriate 
choice to use in an implementation since it is likely to have large coefficients. 
They therefore propose to construct a random vector with small coefficients 
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and check whether (b?-, A) 7^ for all i and j. (Or rather (by, A) 7^ 0, since 
m = in their case.) Only after a fixed number of failed attempts would the 
implementation fall back onto the moment curve. 

Both of these strategies have the disadvantage however that all the terms 
in fll.3p need to be available before the constant term of the first term can 
be computed. This may induce a large memory bottleneck. The authors of 
|DLHTY 04j have therefore also implemented an alternative strategy where 
a random vector with larger coefficients is constructed at the beginning of 
the computation. If the coefficients are large enough, then the probability 
of having constructed an incorrect vector is close to zero. The disadvantage 
of this technique is that the coefficients are larger and that the computation 
has to be redone completely in the unlikely event the vector was incorrect. 

We propose a different strategy which does not require all terms to be 
available, nor does it require the use of large coefficients. We simply repeat- 
edly apply Lemma 12.151 for m! from d — 1 down to m. In each application, 
we can simply use A = 1, which is known to be valid in any case. □ 



We summarize the proof of Proposition 11.111 

Proof of Proposition Given a parametric polytope P C Q n x Q d , apply 
Proposition ^. 41 to obtain the decomposition { Qi }. For each region Qj, apply 
Proposition 12.111 to obtain the corresponding generating function of P s , for 
s G Qi. The result is a collection of polyhedral regions Qi such that, for 
s G int(Qi) n Z n , 

x Pj (s) 

/(P s nZ d ;x) ^ 



(1 - x u ;i)(i - x u ^) ••• (1 - x u ^) 



where Uji G Z d \ {0} and the coordinates of pj : Z n — ► Z d are piecewise 
step-polynomials of degree one. All that remains is to use Lemma 12.151 with 
m = to compute c(s) := f(P s H Z d ; 1) as a step-polynomial in s, valid for 
s G int Qi. □ 
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3 Equivalence of Rational Generating Func- 
tions and Piecewise Step-Polynomials 

In this section, we prove Theorem II .51 that we may convert between rational 
generating function and piecewise step-polynomial representations in poly- 
nomial time. In both directions, we reduce the problem to a set of counting 
problems to which we apply either Proposition 11.111 or Proposition 11.121 We 



first prove a special case of the first half of Theorem 11.51 as a corollary of 
Proposition 11.111 

Corollary 3.1. Fix d. There is a polynomial time algorithm which, given 
a E Q, p G Z n , ai, a 2 , . . . , a^ e Z™\{0} and given 1 6 Z" such that (1, a;) ^ 
for all i, computes the piecewise step-polynomial c : Z n — > Q such that 



is convergent on a neighborhood of e . 

Proof. We may assume, without loss of generality, that (1, a«) < for all i. 
Otherwise, if (1, a*) > for some i, we would apply the identity 



It suffices to prove this corollary for a = 1 and p = 0, because if c'(s) is 
a piecewise step-polynomial representation of the generating function g(x), 
then a ■ c'(s — p) is a piecewise step-polynomial representation of ax p g(x)). 
Note that a = 1 and p = mean that c(s) is the vector partition function 
defined in Example 11.81 

We expand /(x) as a product of infinite geometric series, 




1 



— x 



(3.2) 



1 - x 



d 



/(x) = na+x ai +x 2a * + •■•)• 



Then 



d 



/(e 1 ) = J](l + e ^)+e 2 ^) + -..), 



i=i 
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and this expansion is convergent on a neighborhood of e 1 , since (1, a,) < 0. 
We see that we are looking to compute the function 

c(s) = #{A = (Ai, A 2 , . . . , A d ) G Zi | s = Aiai + A 2 a 2 H h A d a d }. 

Let P be the parametric polytope 



P={(s,A)G 



A > and s = Aiai + ■ — h Xd&d}- 



Then 



c (s) = #{A G Z d | (s,A) G P} 



which can be computed as a piecewise step-polynomial using Proposition !!.!!! 
The proof follows. 

Example 3.3. Consider the function 



□ 



/(x) 



1 



(1 - x( 1 > 1 ))(l - x( 2 - 1 ))(l - x( 1 '°))(l - xCW) ' 
which is the generating function of the vector partition function 

'1 2 1 N 



c(s) =# AGN 



110 1 



(3.4) 



the solution set P s 



A G 



A > and 



is a two 



This is the same as the example from [Be c04i Section 4]. For a given s G Z n , 

'1 2 1 N 
I 1 l y 

dimensional polytope in Q 4 , so it is helpful to convert it to a full-dimensional 
polytope in Q 2 (without changing the number of integer points). To do this, 
extend the transformation matrix from (13 ,4p to 



M 



A 2 1 0\ 
110 1 
10 

\0 1/ 



which is unimodular (that is, it has determinant ±1 and so, as a linear 
transformation, it bijectively maps Z 4 to Z 4 ), and perform the change of 
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coordinates Ah> A' = MA. Then 



c(s) = # A'GZ 4 M _1 A' > and 



12 10 



M -1 A' 



110 1 

# {A' G Z 4 | M _1 A' > and X[ = s h A' 2 = s 2 } 



# (4a;)gz 2 



/ 


-1 


2 


1 -2\ 


( S1 \ 




1 


-1 


-1 1 


S2 










1 


K 


\ 








I J 


W 



> 



This is the enumeration problem that was our running example in the last 
section. 



□ 



We will also need the following lemma. 

Lemma 3.5. Let $(m,d) = (™) + (?) + ••• + (™) ■ Th en m hyperplanes 
in Q d decompose the space into at most $(m, d) polyhedral chambers. Fur- 
thermore, if we fix d, then there is a polynomial time algorithm which, given 
m hyperplanes in Q d , computes the defining inequalities for each of these 
chambers. 

Proof. This lemma is well known, especially the first part, see, for example, 
Section 6.1 of [Mat02] . We prove both parts by induction on m. Certainly the 
statement is true for m—0. Suppose we have a collection of m hyperplanes 
T~(-i, ■ • ■ I'Hmi an d assume that these decompose Q d into at most $(m, d) 
polyhedral chambers whose defining inequalities may be determined in poly- 
nomial time. Let us then add a new hyperplane TC m+ x, which will split some 
of the old chambers in two. The chambers that it splits correspond exactly 
to the chambers that the m hyperplanes TCi D TC m+ i C T~t m +i, for 1 < i < m, 
decompose the (d — l)-dimensional space 7i m +i into. Inductively, there are 
at most <3>(m, d — 1) of these chambers in T~i m +i, an d their descriptions may 
be computed in polynomial time. Therefore, the hyperplanes Hi, ■ ■ ■ ,Tl m +i 
decompose Q m into at most $(m, d) + <3>(m, d — 1) = $(m + 1, d) chambers, 
and we may compute their descriptions in polynomial time. □ 

A generating function in the form (jl.3p is simply the sum of terms like 
those in the statement of Corollary 13 .![ so the first half of Theorem [T75] follows 
from the following lemma. 
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Lemma 3.6. Fix d. There is a polynomial time algorithm which, given piece- 
wise step-polynomials Cj : Z d — > Q, computes c(s) = J2i c i( s ) as a piecewise 
step-polynomial. 

Proof. Suppose Cj(s) are given as piecewise step-polynomials, and let c(s) = 
^Vcj(s). We would like to compute c(s) as a piecewise step-polynomial. 
For each i, let {(ajj,x) < be the collection of linear inequalities that 
define the chambers of the piecewise step-polynomial representation of Cj(s). 
By Lemma I3.5[ we can compute in polynomial time the chambers in Q n 
determined by the collection of all inequalities {(a^x) < These 
will be the chambers in the piecewise step-polynomial representation of c(s). 
Within a particular chamber, each q(s) is defined by 

rii dij 

q(s) = 22 a v II L( a *jfc, s) + b ijk \ , 

3=1 k=l 

where a^- G Q, a ijfe G Q d , and bij k G Q, and so c(s) = ^j Cj(s) is simply a 
sum of such functions. □ 

The first half of Theorem 11.51 is now proved. The second half of Theo- 
rem 11.51 depends on Proposition I1.12[ which we prove now. 

Proof of Proposition \1.12i Given a parametric polytope P C Q n x Q d , apply 
Theorem 4.4 of [BP99J (see Proposition 12. lip directly on P (that is, not 
considering P as a parametric polytope but as a polyhedron in its own right) 
to obtain the rational generating function 

(7(PnZ"+ d ;x,y)= £ x"y* 
(s,t)ePnz™+ d 

in polynomial time. Then the generating function /(x) can be obtained by 
substituting y = 1, i.e., 

/(x)= ^c(s)x s = 9 (FnP +d ;x,l). 

We may perform this substitution in polynomial time using Lemma 12.151 
The result is in the form (jl.3p . □ 
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Given a piecewise step-polynomial c(s), we would like to compute the 
rational generating function /(x) = ^2 seZn c(s)x s . It suffices to prove it for 
functions of the form 



c(s) 



U d j=i l( a j, s ) + b j\ > for s G Q 
0, for s Q, 



where Q is a rational polyhedron, a.,- G Q n , and 6j G Q, because all piecewise 
step-polynomials may be written as linear combinations of functions of this 
form. 

Let P C Q n x Q d be the polyhedron 
P = {( s , t) G Q n x Q d | s G Q and 1 < tj < (a,-, s) + fy, for 1 < j < d). 
Then 

c(s) = #{teZ d |(s,t)GP}, 

and we may compute /(x) = c(s)x s as a rational generating function 
using Proposition 11.121 The second half of Theorem 11.51 follows. 
Finally, we prove Corollary 11.131 

Proof of Corollary \1. left Let 

S = {(s,t) G Z n x Z d | 3u G Z m : (s,t,u) G P} . 

Then we may compute, in polynomial time, the generating function 

/(5;x,y)= ^ xV, 
(s,t)es 

using Theorem 1.7 of [B W03j . Next we compute f(S; x, l) using Lemma [2. 151 
and the c(s) that we desire to compute is the piecewise step-polynomial rep- 
resentation of this generating function. Applying Theorem 11.51 the proof 
follows (since P is bounded, ^ s c ( s ) xS converges everywhere to /(S'jx, l), 
and so any 1 not orthogonal to any of can be used in the application of 
this theorem). □ 
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